Quantum Larmor radiation in conformally flat universe 



Rampei Kimura, Gen Nakamura, Kazuhiro Yamamoto 
Department of Physical Science, Hiroshima University, Higashi- Hiroshima 739-8526, Japan 

We investigate the quantum effect on the Larmor radiation from a moving charge in an expanding 
universe based on the framework of the scalar quantum electrodynamics (SQED). A theoretical 
formula for the radiation energy is derived at the lowest order of the perturbation theory with 
respect to the coupling constant of the SQED. We evaluate the radiation energy on the background 
universe so that the Minkowski spacetime transits to the Milne universe, in which the equation of 
motion for the mode function of the free complex scalar field can be exactly solved in an analytic 

t-H ' way. Then, the result is compared with the WKB approach, in which the equation of motion of the 

t-H , mode function is constructed with the WKB approximation which is valid as long as the Compton 

wavelength is shorter than the Hubble horizon length. This demonstrates that the quantum effect 

O^l ' on the Larmor radiation of the order e 2 h is determined by a non-local integration in time depending 

f— i , on the background expansion. We also compare our result with a recent work by Higuchi and Walker 

[Phys. Rev. D80 105019 (2009)], which investigated the quantum correction to the Larmor radiation 

h— j ' from a charged particle in a non-relativistic motion in a homogeneous electric field. 

in : 
<N ■ 

I. INTRODUCTION 

q-i Recently, the quantum electrodynamics (QED) processes in a strong field background attract interest of theoretical 
researchers as well as experimental researchers [l|-l3| . Especially, recent developments of strong laser technique might 
(3J[)! provide us with a chance to test the Schwinger effect, which causes a pair creation phenomenon in the strong laser- 
'— 1 . electric field [H . On the other hand, some authors claimed that the Unruh effect [3, [f| might be detected by measuring 
the radiation from electron beam in a strong laser background @, H| , but it seems to be necessary for more careful 
investigations for the definite conclusion of the relevance [3] . 
■ An analogy between the radiation process from a moving charge in an expanding universe and the classical Larmor 
Q\ [ radiation in the electro- magnetic theory has been investigated ([13], see also [ill . [l2j]). The motion of a massive 
' charged particle in an expanding background spacetime can be regarded as an accelerated motion because the physical 
momentum of the particle changes as the background spacetime expands or contracts. This leads to the Larmor 
. radiation from the charge in an accelerated motion in the classical electro-magnetic theory. In Ref. [13], it was also 
demonstrated that the quantum effect affects the radiation energy using a model on a background universe with 
specific expansion, where the equation of motion of the mode function is solved in an analytic manner exactly. 
In this paper, we focus our investigation on the quantum effect on the Larmor radiation from a moving charge in an 
J> | expanding universe. Recently, Higuchi and Walker investigated the first order quantum effect of the Larmor radiation 
from a moving charge in a homogeneous electric field background [l3[ • In Ref. [l3[ , the motion of the moving charge 
is limited to the non-relativistic case. A generalization of the work of [HI is presented in Ref. @. In the present 
paper, we apply the method developed in Ref. [3] to the radiation from a moving charge in an expanding universe. 

This paper is organized as follows: In Sec. 2, we derive a general formula for the radiation energy of the quantum 
Larmor process in an expanding universe on the basis of the scalar QED. This formula, which is evaluated at the 
lowest order of the perturbation theory with respect to the coupling constant, is obtained by solving the equation of 
motion of the mode function of the complex scalar field in an expanding universe. In Sec. 3, we consider the complex 
scalar field on the background universe that the Minkowski spacetime transits to the Milne universe, on which the 
equation of motion of the mode function is exactly solved in an analytic manner. Then, the radiation energy is 
evaluated by performing an integration numerically. In Sec. 4, we adopt the WKB approach for the mode function, 
which is constructed in an expanded form with respect to K. This approach yields the formula of the radiation energy 
in an expanded form with respect to h~ as well. The quantum effect of the order of h is determined by the non-local 
nature in time of the background expansion. The result is compared with the work [l3j . Section 5 is devoted to the 
summary and conclusions. In Appendix A, we summarized useful formulas in our derivation of the radiation formula. 
Through out this paper, we set the speed of light c = 1 and the metric convention (+,—,—,—). 
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II. FORMULATION 



First, we derive the formula for the radiation energy from a moving charge in the spatially flat Friedmann- Robertson- 
Walker spacetime, whose line element is given by 



ds 2 = a 2 (rj)[drj 2 — 5ijdx l dx-'\ 1 



(1) 



where 77 is the conformal time and a{r]) is the scale factor. We consider the scalar QED action conformally coupled 
to the curvature 



S = I drjdx^—g 



ieA,. 



1 

44M> 



(2) 



where F^ v is the held strength, /i is the magnetic permeability of vacuum, and £ = 1/6. Introducing the conformally 
rescaled held $, 



we may rewrite the action as 



S = I drjdx 



(D^YD^ -j 

n 



mV^ 1 



(3) 



(4) 



where = t?„ + ieA^/fi. It might be worthwhile to note that the equivalence between Eqs. @ and (U]) will break 
down at one- loop order [ll|. However, this does not affect our result because our investigation is at tree level. 

We adopt the in-in formalism 14] to evaluate the radiation energy, in which the expectation value of the Hciscnbcrg 
operator Q at the time 77 is given by the following form, 



00 N „ n „ V2 

(Q(v)) = J2tn dm--- d m ([if/M, . . . [H^m), Q( v )} ■ - .]> 

J\J- — Q M J — OO J— OO 



(5) 



where Hi denotes the interaction Hamiltonian operator, 



Hi(t]) = -~ I d 3 x^(ry,x)[9 A1 $t( r/;X ) $(?7iX) _ $ t( ?/;X )^ i$ ( ?7iX )] j 



(6) 



and the operators in the right hand side of Eq. ([5]) are those of the interaction picture. The free field of the electro- 
magnetic field quantized in a finite box is 



A=l,2 k 



,A e -ifet, + a ^ e ikv 



(7) 



where denotes the polarization tensor, A denotes the direction of polarization, V is the volume of the box, and 
a£ are the creation and annihilation operator, respectively, which satisfy the following commutation relation, 



In the similar way, the quantized complex scalar field is given by 

T 



(8) 



(9) 



where the mode function f p (r]) obeys the equation of motion, 



d 2 p 2 + m 2 a 2 (ri) 



drj 2 



<P P (V) = Q : 



(10) 
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FIG. 1: Feynman Diagram for the process. 



with the normalization condition ip p tp p ~ (p P <p P — i, and the creation and annihilation operators satisfy the commu- 
tation relation, 



K^U = <W: [ c P' c p'] = W> [OTHER] = 0. 



(11) 



Choosing the Heisenberg operator Q as the number operator a^afc, we calculate the total radiation energy as 
follows, 



E= d 3 khkJ2(ak a t) =~ d 3 khkJ2 h ~ 2 / *a / *i ( [Hi(h), [F f (*a), a^'aj^ 

«/ \ «/ \ </ — CO — OO 

where we only considered the contribution up to the order of e 2 . We set the initial state as, 

|*> = |0)„®|1)6®|0) 5) 



(12) 



(13) 



where the subscript a means the photon number basis, b is the charged particle number basis, and b is the anti-particle 
number basis. 

We consider the lowest order contribution of the process so that one photon is emitted from the charged particle, 
as shown in Fig. [TJ From Eq. (|T2|). we have 



d 3 k 
e J (2k) 3 



(^ P/(?7)Vpi(77) " ViM^PiW) 
dT, ~/2k \~ (fipf ^ + ^Pf (^) X^ Pl ^ 



(14) 



where Pf = Pi — Kk, Kk is the photon momentum, and eo is the permittivity of vacuum, which is related to /zq by 
eoMo = l/c 2 = 1- Using the equation of motion (|10p. Eq. (|14p leads to 



E 



d 3 k 
eoJ W 



2fc ft 
e tkri gi + p/ 



where k = k/|k|. Furthermore, Eq. (fl5|) reduces to (cf. [To| ) 



E = 



d 3 k 4(p? - (pi • k) 2 ) 



2e J (2tt) 



(15) 



(16) 



Note that this is the energy in the conformally rescaled spacetime, which is not the physical energy. We can read the 
physical radiation energy as £ = E/a. 



III. SOLVABLE MODEL IN MINKOWSKI-MILNE UNIVERSE 



In this section, we consider the universe whose scale factor is written in the form, 

a(r)) = y/l + e 2 w, 



(17) 
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where p is the parameter to describe the expansion rate. This gives the Minkowski spacetime at r\ — > — oo, and the 
Milne universe at r\ — > +00. The equation of motion of the mode function (|10[) is 



d 2 p 2 +m 2 (l + e 2pr i) 



_d V 2 h 2 
The solution is given by using the Bessel function, 



Mv) = 0. (18) 



where 



/ p 2 + m 2 _ m 

P 2 h 2 ' m =p~h- 

The radiation energy is given by the formula (|16p , then we focus on 



/ 7T / 7T fh~ ik IP r°° 

Y 2psinh< f V 2psinh7r& p 7 f 



where we defined 



p 2 +m 2 pj + m 2 — 2hkpi cos 9 + h 2 k 2 + m 2 

^ = VZ^' e/ = VZ^V 7* ■ (2) 



We use the mathematical formula I If 



Mj ( z)J^ ( z) rm _(.._ K)/2+ l)vm + f< _ K)/2 + i )r( z(_0 + 6 - «)/2 + 1) : 
where we defined k = k/p. Then, the radiation energy is 



(22) 



2 9 

E = Fl 



2ttc p' 



J dkk 2 J d cos 6»(1 -cos 2 6) 



sinh7r£/ sinh7r£j sinh7rft (£/ — + k) shm7r(£/ — £j + K )/2 
sinh7r(£/ — k)/2 sinh7r(£/ + & + k)/2 sinli7r(£/ — & — k)/2 
• 6 £/ + 6 + « £/-&-« 



(23) 



The classical radiation energy corresponding to the Larmor radiation can be evaluated using Eq. (|4"3"f (see also 
Eq. (36) of Ref. pjj). In the case of the background universe with Eq. (TiTl) . we have 



247reom 2 

Combining (l23l) and (p4|) . we have 



= (24) 



E c i 



12 f ^ dux 2 J dcos8(l - cos 2 



sinli7r£/ sinh7r^ sinh7TK (£/ — £j + k) sinh7r(£/ — 6 + K )/2 
^ sinhTr^/ + & - k)/2 sinh + £ + /c)/2 sinliTr^/ - & - k)/2 

+ + & + « 

We performed the integration of Eq. ()25[) numerically. Figure [2] plots E/E& as a function of m/pTi with fixing 
Pi/m = 0.1 (black solid curve), Pi/m = 1 (red dashed curve), and Pi/m = 2 (blue dotted curve). In the limit 
m/ph 1, one can see that E/E c \ approaches to 1. On the other hand, in the limit m/ph <C 1, E/E c \ approaches to 
0. Complimentary to this figure, Figure [3] shows the contour of E/E c \ on the plane Pi/m and mj 'ph. For the region 
m/ph 3> 1, E/E c \ approaches to 1, while for the region m/ph <C 1, E/E c \ approaches to 0. These features are the 
same as those in the previous paper (loj . though a different background universe was considered there. 
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FIG. 2: The ratio of the total radiation energy to the classical WKB formula, E/E c \, as function of m/ph with fixing Pi/m = 0.1 
(black solid curve), Pi/m = 1 (red dashed curve), and Pi/m — 2 (blue dotted curve). 
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FIG. 3: The deviation from the classical WKB formula. Contour of E/E c \ on the Pi/m and m/ph plane. 



IV. WKB APPROACH 

In the previous section, we investigated the radiation energy in the case when the mode function is exactly solved 
in the analytic manner. From the result of the exactly solvable model in the previous section as well as the previous 
paper [10| . the quantum correction reduces the total radiation energy compared with the classical Larmor formula. 
In this section, we consider the origin of the quantum effect in the Larmor radiation. To this end, we adopt the WKB 
approach. Our computation is a generalization of the work by Higuchi and Walker [l3| . in which the quantum effect 
of the Larmor radiation on an electric- field background is investigated (cf. Q). 

The mode function with the WKB approximation is 



1 



exp 



% \ n p {rj')dr}' 



(26) 
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with 



\/p 2 + m 2 a 2 (r)) 



(27) 



Substituting the WKB solution (|26|) into Eq. (fl"6f . we have 



E 



d 3 k 



e 2 1 



2e 2 2 y (2tt) 



k • (2p j; - /ik) 
drj — — j== — 7= exp 



ifc?y + i / (f2 P/ — f2 p Jd?/ 



2p t - Kk 
drj == — -= exp 



ihj + i / (f2 p/ — Q Pi )drj' 



with 



^p* fa) = V P; f + ™ 2 a 2 fa) A, 



fi p/ fa) = V(P* - Kk) 2 +m 2 a 2 (r 1 )/K 

We derive the expression for the radiation energy in the form expanded as a power series with respect to K. 
the approximation, 



(28) 

(29) 
(30) 
Using 



i 



k • p t 



(k ■ Pl ) 2 



y/p 2 +m 2 a 2 ( V ) 2 \ y /p 2 + m 2 a 2 (7 1 ) ^ p 2 + m 2 a 2 ^ 



1 



k • Pi 



hy/n^y/n^ y/p 2 + m 2 a 2 ( V ) r ' 2p 2 +m 2 a 2 ( V ) 
we have the expression of the radiation energy up to the order of h, 

2e 
1 



2e J (2tt) 3 

Tik k • pi 



d£/V d£'/ d£ df 
k • pi 



Jk£-ik£' 



2 p 2 +m 2 a 2 (r]) 2 p 2 + m 2 a 2 (?/) 
in /"V fc 2 fc 2 (k- Pl ) 2 



2 J n , V ^p 2 +m 2 a 2 {i 1 ") ,/p 2 +m 2 a 2 (r ] "y 
where we introduced the new variable £ instead of r\ 

rv k • p, 



drj" + 0(h 2 ) 



t = v- 



VPi + rri 2 a 2 {il') 



drf. 



Furthermore, we introduce the notation, 



dr 
dr\ 



Pi, 



\/p 2 + m 2 a 2 (ri). 



Then, Eq. ([51]) is rephrased as 





e 2 


/■ d 3 k 


2e 


1 (2^) 3 








X 


1 + 


»( t . 









dx dr ^ dx' dr' 
drj drj drj 1 drj' 



ink 2 r (s) 



dr" 1 



(f) 



dx" 



(31) 

(32) 

(33) 
(34) 



(35) 
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We replace k in Eq. (|35|) with the partial derivative with respect to £ or £' operating on e lk ^ , which is the 
equivalent technique adopted in Ref. [13l |. Partial integrations lead to 



E (o) 



2e (2 

„2 



/^/«/^>(( t .£)(^~)_£~), 



(36) 



2e (27r) 3 
\ ih ( d 



r/k 



dk d£ d£e 



d_d_ 
4 \d£ d(' J d£d£' 

ih d 2 d 2 r 



/ - dx \ / - dx' \ dx dx' \ I - c(x ctr - dx' o(t' 

V "d£/V 'df 7 / ~ ~dl "de ) [ "ch~ch + "dn'dn' 



dx' 



dx dx.' 



dx dr 
drj drj 



dx' dr' 
drj' drj 1 



2 d£ 2 d£' 2 



dx' 



dx dx' 



1- k- 



dx" 

d 7 / 



(37) 



where i?( ) and E^ are the contribution of the order of h° and h, respectively. The expression of the order of h°, 
Eq. dSH]), gives the classical radiation energy formula 

E3 EL(cf- EMI), whi le the expression of the order of h, 
Eq. (|37l) , gives the first order quantum effect in the radiation [9|, [l3| . We choose the coordinate z to be parallel to pi , 
and the polar angle 9 to be the angle between the z-axis and photon momentum vector k. Then, we write 



dz 
dry' 



— = (1 — ?jcos6>), 
drj 

d 2 z _ v 

~d~e ~ 



(1 — V cos ( 



with which we have 



E (0) 



/ d cos 6>(1 -cos 2 6>) / d£ 
87re 7_i J 



d 2 z 

~d~e 



87re J_ 

from Eq. (|36p . Here, the 'dof denotes the differentiation with respect to rj. Integration with respect to 9 yields 

E (o) 



67T60 



drj 



(l-v 2 ) 3 ' 



With the use of the relation v — Pi/ \/p 2 + m 2 d 2 (rj), Eq. (|4"2l reduces to 



=,2„2 



E (0) = _±Ei 
67reo?n 



dr]—. 



(38) 
(39) 
(40) 

(41) 
(42) 
(43) 



This is the classical radiation formula in an expanding universe which can be regarded as the Larmor radiation fic| . 

Next, we consider the contribution of the order of h, Eq. (|37p . Let us first consider the limit of the non-relativistic 
motion, u«l, we can write the radiation energy in terms of the velocity and its time derivatives, and find the leading 
contribution (see Eq. (|A3|) in Appendix A), 



3 /A 2 [ dcos9 [ dr] ( ' dq'—^—j 
o(2tt)^ J J J rj-r)' 

((vv'v' — v'vv) + i>v'(3(vv' — v'v) — (vv — u'u')) 008 6*1 



1 — cos 2 



(44) 



where the 'dash' denotes the function of rf , i.e., v' = v(r)'). Note that the latter terms of the right- hand- side of 
Eq. (|44l) . which is in proportion to cos 9, give no contribution when integrated over 9. In the non-relativistic case, we 
have 
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This gives a rough estimation of the ratio of the first term (vv'v' — v'vv) in the right-hand-side of Eq. (l44l) to the 
latter terms in proportion to cos#. Therefore, the first term is larger than the latter terms in the non-relativistic case, 
ma/pi ^> 1. 

Next, we consider the relativistic limit, v ~ 1. In an expanding universe, the condition of the relativistic motion is 
guaranteed as long as ma(rj) £ p,. Hence, the particle motion becomes non-relativistic as the universe expands with 
(|17[). even though the motion is relativistic initially. In the model with (fT7|) . the particle motion is relativistic for 



Pi 

On the other hand, the ratio of v to v 2 is written as 



< 1. 



(46) 



V 




da + a 2 


V 2 


m 2 


(da) 2 



Pj_ -2p V 
7 



(47) 



This means that \v 2 \ <C \v\ while the particle motion is relativistic. In this case, using expression (IA5[1 in Appendix 
A, we have 



3e 2 h 
4e (27r) 2 
~ 1 — cos 2 



dcosO drj drj 



1 



pi 



(1 — v cos 6*) 3 (1 — v' cos ( 



(1 - v cos 6) 2 (1 -v' cos 0) 3 



(48) 



This expression is derived assuming the particle motion is relativistic all the time. Then, it is not justified when the 
particle motion changes from relativistic to non-relativistic as the universe expands, as for the model with (1171) in 
Sec. 3. 



V. DISCUSSIONS 



Integrating the right-hand-side of Eq. (|44j) with respect to 0, we obtain the quantum correction to the Larmor 
radiation of the order of h, 

= — 2 / d V / drf — . (49) 

47r 2 e Pi J J T]-r)' 

This result should be compared with the result by Higuchi and Walker (131 ] , which presented the quantum correction 
to the Larmor radiation of the order of % from a charged particle in an accelerated motion on a homogeneous electric 
field background, 

,-,fi1 g 2 T% f , f ,,vv' — v'v 

E ™ = i^J dt J dt ^ir- (50) 

There is a difference between (|49| and (|50|) . i.e., the numerical coefficient and the dependence on the velocity and its 
derivatives. However, the quantum effect is non-local in time in both cases. 
In the case when the velocity v and its derivatives are approximated as 

Pi Pi /K-i \ 

V = ~ T— T, (51) 

^p 2 + m 2 a 2 (ri) ma{r)) 

v ~ 2 ' ( 52 ) 

Pi da — 2d 2 



m 

Eq. g9l) yields 



(53) 



E(1) = e^pl I r ^_ (aa- 2^-^-2^ 
47T z eom m z J J rj — rj \ a J a' J 



where a = a(rf) and a' = a{rj'). 
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FIG. 4: Comparison of the deviation of the classical radiation from the total radiation. The solid (red) curve is —E^'/E c \ as 
a function of m/ph, while the other curves are —{E — E c \)/E c \ with fixing pi/m — 0.01 (black dashed curve), pi/m = 1 (blue 
dotted curve), and Pi/m = 2 (green dash-dotted curve), where E^ is computed using Eq. (|54[) . 



Figure [4] shows the deviation of the total Larmor radiation from the classical Larmor radiation from a moving 
charge in the universe considered in Sec. 3, as a function of m/ph. The curves compare the case when the radiation 
is computed with the exact mode function of Sec. 3 and the case when the radiation is computed with the WKB 
approaches up to the order of h. The (black) dashed curve, the (blue) dotted curve, and the (green) dash-dotted 
curve are the results with the exact mode function with Pi/m = 0.01, 1, and 2, respectively. The solid (red) curve 
is from the WKB mode function, Eq. (|54j) . For the region with Pi/m <C 1 and m/ph 3> 1, the quantum correction 
can be well explained by the WKB approximate formula, Eq. ([Ml) . The condition Pi/m -C 1 comes from the fact 
that Eq. (|54p is derived in the non-relativistic approximation. The deviation of the results with pi /m = 1 and 2 from 
Eq. ([54| comes from the fact that the particle motion is relativistic initially. However, the deviation is not significant, 
which may indicate that the fraction of the radiation from the particle in the non-relativistic regime is not small. 

Finally, we show that the condition m/ph 3> 1 is the necessary condition for the WKB approximation in an explicit 
manner. The WKB approximation holds under the condition |19j |. 



1 

2^2 



19l 

2ft 2 



< 1. 



(55) 



In the universe with the scale factor (fl7)) . the condition (|55|) is written as 



hp 



(a 2 -l)(a 2 -5-4p 2 /m 2 ) 



4(a 2 + pf/m 2 ) 3 



< 1. 



(56) 



Hence, the use of the WKB approximation is guaranteed when m/ph 3> 1 is satisfied. This condition is interpreted 
as follows (lo| . In the universe with p7|) . the ratio of the Compton wavelength to the Hubble horizon is 



hp a 2 — 1 



(57) 



When the condition m/ph 3> 1 is satisfied, the ratio (1571) is much less than unity. Thus, the condition for the WKB 
approximation can be regarded as the condition that the Compton wavelength is much shorter than the Hubble 
horizon. 



VI. CONCLUSIONS 



In this paper, we investigated the quantum effect of the Larmor radiation from a moving charge in an expanding 
universe based on the framework of the perturbative approach of the scalar QED theory. We first derived the general 
formula for the quantum radiation energy, which is expressed by the mode function of the free complex scalar field. 
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Assuming the background universe that the Minkowski spacetime transits to the Milne universe, in which the mode 
function is exactly solved in the analytic manner, we demonstrated how the total radiation energy deviates from the 
classical Larmor radiation formula. This gives us a hint to understand the quantum effect on the Larmor radiation. 

In the latter part of the paper, we investigated the quantum effect of the Larmor radiation with the use of the 
WKB approximation in constructing the mode function. In the WKB approach, which is valid as long as the Compton 
wavelength is shorter than the Hubble horizon length, we investigated the quantum effect of the Larmor radiation 
by evaluating the contribution of the order of fi of the radiation energy with the approximate mode function. The 
radiation of the order of fi well explains the deviation of the total radiation with the exact approach from the classical 
Larmor formula in Sec. 3. The radiation of the order of fi is determined by the non-local integration in time depending 
the behavior of the expansion of the background universe, as in the case of a homogeneous electric field background 
0, [l3j]. However, the first order quantum correction to the Larmor radiation in the conformally flat universe is not 
the same as that on the homogeneous electric field background. 
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In this Appendix, we summarize useful formulas in deriving the radiation formula in the present paper. Using 
Eqs. ([33]), flM]), and (|38])-(|40]), we have 



Appendix A: useful formulas 



dt? 



(1 - V COS0) 4 



V 



+ 



(1 — VCOSd) 5 ' 



3v 2 cos 8 




de 




+ 



Pi(l — v cos 8) 



2vv 



2 - 




11 



where the dot denotes the differentiation with respect to the conformal time rj. Using these relations and k • dx/dry 
cos Odz/drj = vcosO, it is straightforward to derive 



d d \ d d 
di ~ ~d£) d£ df 7 



/ - dx \ / - dx' \ dx dx' \ ( - dx dT - dx' dT' 
V 'd?Jl 'dFv ~ df "dT 7 \ 'd^ + 'drfdrf 



dx 



(cos 2 6 — 1) cos ( 



+ 



dx dx' 
df ' df 

3v 2 cos 



dx dr 



dx' dr' 



d?7 dr\ 



drf dv[ 



(1-v cos 6>) 4 (1 - d cos Of J \ {l-v' cos 6>) 3 
w' 3-i/ 2 cos 6> 



(1 — V COS 

!> 



1 — v cos 6* 

£> 

(1 — i>cos( 
v 

(1 — i>cos( 
v 

(1 — ucos( 



(1 — vcosO 
v' 

(1 — v' COS 



(l-v'cos6>) 4 (l-w'cos6») 5 
2vv 2v'v' 



Pi(l — vcos9) pi(l — v'cos( 



(1 - t/cosfl) 3 
v' 

1 — v' cos 8 
3v 2 cos 9 
(1 — v cos 9] 

v' 

(1 — v' cos 9 



2v 2 



(1 — v cos 9) 3 

2v' 2 
(1 — t>' cos 9) 
v' 

(1 — v' cost 
3v' 2 cos 9 



2vv \ 1 
(1 — vcos9) 2 J pi 

2v'v' \ 1 
(1 -v'cos9) 2 

2v'v' 
Pi(l — v' cost 

2vv 



9) 2 )p, 



(1 — v 1 cos 9) 5 ) pi(l — v cos 



(Al) 



and 



d 2 d 2 



di 2 d? 2 

= (cos 2 9-1) 

- 2 



dx'^ 

df 



fd 3 zd 3 z' 



\d£ 3 d£' ; 



dx dx' 

d£ 



5 d£" 



+ 



3w 2 cos( 



d£" 



(»?') 



1 — cos 



1 



dx" 
d^ 77 



dz'V N 
d?7"y 

+ i>' cos6») 



(1 - vcostf) 4 (1 - ucosfl) 5 / V(l - v'cos9) 
v' 3v' 2 cos 9 



(1 — v cost 
v 



(l-v'cos9) 4 (1 — v'cost 
3v 2 cos 9 



(1 — wcos6*) 4 (1 — vcos( 



Pi 

v(l + vcost 
Pi 

i'(l + 2v'cos9) 



v 



1 — v cos t 



1 — v' cos 9 J pi (1 — v' cos t 
3i>' 2 cos9 \ v(l + 2vcos( 



(1 — u'cos^) 4 (1 — v' cos9) 5 J Pi(l — vcos 



In the limit of the non-relativistic motion, «<1, the above formulas give 



ih 

T' 



d d \ d d 
dff ~ df J d^de 

d 2 d 2 



dx^ 



d£ 2 d£> 2 



dx dx' 



dx dr 
drj drj 



dx' dT' 
drj' drj' 



{ k -dc){ k -de)-dz-d?j J ew) d t dr y 1 [ k ■ W 



3ih 
~2~ 



1 — cos 



Pi 



{vv'v' — v vv + vv (3(w — it;) — (vv — v v )) cos 9} 



Here, we omitted the contribution from the term, 

^0-1)^— — 



4 „„dr" ( 9n /dz"\ 2 
d£'—\l-cos 2 9[ 



d£" 



\drf' 



(A2) 



(A3) 



(A4) 



because its leading contribution to the radiation energy in Eq. (|37[) is zero when v is zero at the boundary r\ — , ±oo. 
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In the limit of the relativistic motion, v 



1, we have 




d_ _d_ 



dx 



k- 



dx' 



d£ ) \ d£> 



/ d 2 
l d£ 2 d£< 2 



ih 



.,2 , 



~d? 



dx dx 



dx dx' 
dt; ' 



dt di 



k 



At" 



dx dr 
drj drj 



k- 



'W) 



d£," 



dx' dr' 
drj' drj' 

.11 \ 2 



dr]" 



4 pi(l — v cos0) 2 (l — v' cos( 



v v 



(1 — v cos 9) 2 {1 — v 1 cos 9) (1 — v cos 9)(1 — w' cos ( 



(A5) 



Here, again we omitted the contribution from the term (|A4[) because of the similar reason to that in the non-relativistic 
limit, as long as the particle motion is relativistic all the time. However, the omission is not justified when the particle 
motion changes from relativistic to non-relativistic as the universe expands. 
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